The Markov approximation for the atomic output coupler 
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i The regions of validity of the Markov approximation for the coupling of atoms out of an atomic 

. trap are determined. We consider radio-frequency output coupling in the presence of gravity and 

collisional repulsion, and Raman output coupling. The Markov approximation is crucial in most 
theoretical descriptions of an atom laser that assume a continuous process of output coupling from a 
trapped Bose-Einstein condensate. In this regime many techniques proved to be useful for modeling 
the optical laser, such as master equations, can be used to describe the dynamics of the damping of 
^ C"| the condensate mode undergoing output coupling. 

o 

£ 

I. INTRODUCTION 

The recent observation of Bose-Einstein condensation (BEC) in atomic traps has attracted widespread atten- 
tion. One of the most important perspectives of this experimental achievement is the possibility of producing the 
I | matter wave analog of a laser, i.e., a high flux source of coherent atoms. In its simplest form, such an atom laser can 
be built by adding a suitable output coupling mechanism [[||4|] to present Bose-Einstein condensation experiments. 
Several attempts have already been made to develop a theoretical description of an atom laser |3|-|l7| that combines 
elements of kinetic theory and laser theory. One of the problems that exists in applying optical laser concepts to 
the Bose-Einstein output coupler situation is the question of the validity of the Markov approximation The 
Markov approximation is an extremely powerful tool to describe the coupling of a system of trapped particles to an 
environment. It allows one to think of the coupling in the following terms: at any one time, a particle is either in the 
trapped system or has been coupled out. Quantum mechanically, this means that the existence of any superposition 
of these two possibilities is neglected. Such an assumption is valid if the superposition decays on a time scale much 
faster than changes occur in the state of the trapped system. Classically, the Markov approximation implies that 
an atom that has been coupled out will have no chance of being brought back into the trap again. In the case of a 
7-H ' noninteracting gas and ignoring gravity, atoms that are coupled out leave the spatial region of the trapped atoms due 
to the relatively slow quantum spreading of their wave packet. As a consequence, the superposition decays slowly 
and it is likely that some atoms will be coupled back into the trap instead of leaving the system irretrievably. Such 
a behavior leads to a strongly non-Markovian dynamics of the condensate mode 0-0]. However, in the presence of 
accelerating potentials atoms may be removed from the region of coupling at a much faster rate and thereby allow 
the use of the Markov approximation. 
. In general, there are two distinct operating regimes for an output coupler. The first regime is the strong-coupling 
regime, where portions of a condensate are coupled out at such a rate that the output-coupled atoms do not have 
time to propagate while the coupling is in progress. This regime was realized in the first experimental demonstration 
O ■ of an output coupler j^] , where a rf-pulse was used to couple out a large fraction of the trapped atoms within a time 
interval of the order of fis. When the coupling time was long compared to the rate at which atoms can be coupled 
back into the trap, strongly non-Markovian behavior, such as Rabi oscillations, was observed. In this regime strong 
collisional interactions dominated the dynamics of the output-coupled atoms. The Gross-Pitaevskii equation (GPE) 
has been found useful in numerically modeling such a situation p-|ll|. However, in this paper we are interested in 
the opposite limit of weak semi-continuous output coupling, where the untrapped atom beam is of low enough spatial 
density that we can neglect the effect of collisions on its behavior. Output coupling in this regime has not yet been 
demonstrated, although it is considered experimentally viable j3|. 

Previous theoretical treatments have implicitly made use of the Markov approximation by assuming a Lindblad 
master equation (see Ref. |l8|]) for the condensate mode JI^-pj|. A different, though essentially equivalent, approach 
treats the loss of condensate atoms by adding a damping term to the time-dependent Gross-Pitaevskii equation 
p6| , pT| . In the present work we explore the regimes where the Markov approximation is valid. In Sec. [ij we present 
the Hamiltonian of the total system and introduce the concepts of a memory function and a memory time. In Sec. 
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Ill we apply the ge neral ideas to the specific case of the radio-frequency output coupler. The results obtained are 
summarized in Sec. [II D . In Sec. IV we also address the case of a Raman output coupler. 



II. MATHEMATICAL FORMALISM 



The Hamiltonian of the total system can be written as the sum of three parts 

H = Ht + Hjj + Hi, 



(1) 



where Ht is the Hamiltonian of the trapped particles, Hjj is the Hamiltonian of the untrapped particles and Hi the 
interaction between the two. The Hamiltonian of the trapped atoms is of the form 



H T = / dx V>t( x ) 



2m 



^t(x) + ^ / dx V4( x )V't( x )^t(x)V't(x), 



(2) 



where ?At(x) is the field operator for the trapped atoms and Vr(x) is the trap potential, m is the mass of the atoms 
and U — 4Trh 2 a/m is the coupling constant for a local collision process, where a is the scattering length of a trapped- 
trapped collision. The untrapped atoms are assumed to be of low enough spatial density that we can ignore collisions 
between the untrapped atoms and can write 



Hjj = I dx^(x) 



2m 



(3) 



where V>t/(x) is the field operator of the untrapped atoms and V{} (x) is the effective non-confining potential experi- 
enced by the untrapped atoms. The form of this effective potential will be given in a later section. Both fields satisfy 
Bose commutation relations, [^(x), = 5(x — x'). We consider situations where the output coupling is linear 
and conserves particle number, 

Hi = ihy^ J dx [$(x,t)$,(x)$r(x) - ff*(x,t)^(x)^(x)] , (4) 

where the coupling constant, g(x, <), is normalized so that J dx|g(x, t)\ 2 = 1, and we assume that the strength 
of the coupling, yfiy, is time independent. The time dependence of g(x, t) is then simply an oscillatory phase, 
g(x,t) = g(pc)e~ ivt . The above interaction Hamiltonian can describe radio-frequency output coupling ||, where a 
radio wave induces a transition from a hyperfine level that is trapped in the magnetic trap to one that is untrapped 
or anti-trapped. It can also describe Raman output coupling where two laser beams cause an atom in the trap to 
make a two-photon transition to an untrapped state. In this case the atom experiences a momentum kick equal to 
the difference in the momentum of the photons involved in the transition. 

In this work we confine our interest to a single energy mode of the trapped system. We will concentrate on the case 
when the trapped mode of interest is the condensate mode. However, the method is equally applicable to an excited 
mode of the trap. The trap mode operator of interest is defined in terms of its spatial mode function, u a (x), by 



dx u* a (x)tp T {x). 



(5) 



An oscillation frequency, n, will be associated with this mode. If the mode is the condensate mode then u Q (x) is 
the solution to the time-independent GPE and fi is the chemical potential. The mode operator a is coupled to the 
untrapped field throughout an effective coupling region given by 



K(x,t) = g(x.,t)ua(x). 



(6) 



A. Damping of the trapped mode 



By substituting the formal solution of the Hciscnbcrg equation for the untrapped field into the Heisenberg equation 
for the trapped mode a, we obtain the Langevin equation of motion [19], 



2 



^L = - l -[a{t),H T ]- 1 daf m (t-s)a(a)-Syt(t)- (?) 

The driving field, is the contribution from the free propagation of the initial untrapped field, considered here to 
be in a vacuum state, 

£(t)= f dx«*(x,t)^°(x,t), (8) 



where the dynamics of the free untrapped field operator ipy are determined by the Hamiltonian H\j alone. The 
presence of the driving field is necessary to preserve the bosonic commutation relations of the mode operator a. The 
damping term (the second term on the right-hand-side) represents a loss of particles from the trapped mode into the 
untrapped field and makes explicit the dependence of the trapped mode on its past behavior via the so called memory 
function, / m (t — t'). The correlation between the driving field and itself at an earlier time determines the memory 
function via the commutation relation, 

/m (t_t') = £t (9 ) 

This relationship between the driving field and the memory function is an example of the quantum fluctuation- 
dissipation relation and leads to a description of a damping process consistent with both quantum and statistical 
mechanics. The damping process can then be interpreted as the coupling of discrete atoms out of the trap at random 
times. 

In general, the Langevin equation (]?]) will contain a term representing a free oscillation and other terms due to 
collisions between the trapped modes. We are interested in output coupling and the collisional behavior of the trapped 
atoms is not explicitly modeled in this work. Instead, we define a rate T to account for these other processes without 
considering them explicitly. This rate may be calculated in work concentrating on the trapped atoms, such as [ p0| . 

B. A finite memory time 

For a dissipative system we expect the system behavior at time s as s — > — oo to become less and less important 
in determining the present behavior of the system. To make this more concrete we define a memory time, T m , as the 
time after which we can neglect the effect of the previous behavior of the system on the evolution in the present. In 
terms of the Langevin equation, a memory time exists for the system if, at some finite time in the past, T m , we can 
make the approximation 



dsf m (t - s)a(s) ~ / dsf m (t - s)a(s), (10) 

> Jt-T m 

for all t. If this condition is satisfied then we will call T m the memory time of the system. Note that if a memory time 
cannot be defined for the system then the separation of the total system into a localized system interacting with an 
environment becomes inappropriate as there will be no time at which one can say that a particle has left the localized 
system and entered the environment. 

In general, the memory time as given by Eq. (|Io| ) depends on the nature of a(s) and so no general statements can be 
made concerning this condition without detailed knowledge of the behavior of the trapped atoms. However, from the 
motion of the untrapped atoms one can determine cases where condition Eq. ( |l0|) can be satisfied without assuming 
too much about the behavior of the trapped system. This is clearer if we write the memory function in terms of the 
single particle Green's function as, 

fm(t-t') = J dxdx' K*(x,t)K(x',t')G(x,t;x',t'), (11) 

where 

G(x,t;x',i') = [^(x,i),^(x',t')], (12) 

= ({0}|^(x,<)V^(x',t')|{0}), (13) 

is the single particle Green's function for the free atoms for t > t' and k(x, t) is the effective interaction region given 
by Eq. (|^) . The memory function can therefore be interpreted as the overlap between an atom, with an initial wave 
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packet of the shape of the interaction region, «(x', t'), with a wave packet, k*(x, t) after it has propagated for a time 
t — if. There are two distinct ways in which this overlap could become smaller with increasing time. Firstly, if an atom 
is leaving the interaction region due to an accelerating potential or quantum mechanical spreading of its wave packet 
then the overlap between the atoms wave packet and its original wave packet will decrease in time. The time at which 
we can neglect this overlap will determine the memory time. Secondly, if an atom is accelerating then it will gain 
kinetic energy which will cause the wave packet of the atom to oscillate. The overlap between this wave packet and 
a stationary one will then also oscillate. This oscillation will average to zero when integrated over time scales much 
longer than the period of oscillation. In our case, we are interested in integrating over time scales short compared to 
the coupling time scale, 7 -1 . The memory time can be defined as the time when the oscillation is much faster than 
7, as the overlap will average to zero after this time. Another, equivalent, way of thinking about this is that only a 
range of frequencies will be close to resonance with the coupling. Particles with energies far from resonance will not 
be coupled back into the trap. 

Following from these considerations we can make a more practical definition of a memory time in terms of the 
motion of the untrapped atoms alone by 



dsf m (t - s) 



> 



t-T m 



dsf m {t - s) 



(14) 



independent of a(t). Note that we have made the replacement f m (t — t') — ► f m (t — t')e~ ip,i to take account of the 
oscillation of the mode ait), as this may cancel the oscillation of the memory function itself. This is the only aspect 
of the behavior of ait) that we will consider in determining a memory time. 

In summary, there is a certain region of phase space in the untrapped field where particles can be coupled back into 
the trap. The time taken for a particle to leave this region of phase space determines the memory time. The memory- 
time can then be interpreted as the time interval after which we can safely assume that a particle has irretrievable 
left the trap. 

If a change in a(t) occurs during the memory time then it is necessary to consider the coupling out of an additional 
particle before the first particle has either been coupled back into the trap or left irretrievably; this is the non- 
Markovian regime. Strictly speaking we should distinguish between the free evolution of a(t) and its evolution due to 
the output coupling. Strong coupling, jT m > 1, leads to the second-order effects mentioned above and this is a serious 
breakdown of the Markov approximation. The neglect of these effects for weak coupling is often referred to as the 
Born approximation. If on top of this, the free system evolution (except for an oscillating phase), is on a time-scale 
much slower than the memory-time, TT m -C 1, then a(s) can be taken to the front of the integral in the damping 
term of the Langevin equation and the integral over the memory function can be done. The equation for a(t) will 
then be local in time; this is referred to as the Markov approximation for the damping. 

In our case, the Markov approximation can be made for the damping if the memory-time, T m , is much less than 
the time-scale of both the evolution of the trapped mode, T _1 and the damping, 7" 1 . The operator a(s) can then be 
replaced by its value at t so that Eq. (Q) becomes, 



da(t) 
dt 



where 



[a(t),H T ] - iAu;a(t) - j'a(t) - ^(t), 



7' = 7^ / dsf m {t-s) 



Ace; = 7 3 jy dsf m (t-s)\ 



(15) 

(16) 
(17) 



where Aui is a frequency shift due to the coupling |18| . Equation (|15|) is equivalent to the master equation for the 
reduced density matrix of the system Jl^] . 

For the condensate mode the oscillation frequency and the decay rate (determined by the spatial mode function) will 
vary at a rate 7 with the number of atoms in the mode. If the collective excitations of the trapped atoms (caused by 
the change in population) decay much more rapidly than the loss rate then the chemical potential and the condensate 
mode function determined from the time-independent GPE will be valid on time scales much shorter than those of the 
loss. Steck et al. |l6| ] have demonstrated that this procedure is valid by numerically simulating the evolution of the 
full coupled GPE for the trapped and untrapped fields. In the present case Eq. §LQ) not only describes the damping 
of the condensate number, it also determines the approximate evolution of the total quantum state of the condensate 
mode (assuming a slow phase diffusion pl|-p3|]). 
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It is likely that a continuous-wave atom laser will also have some form of pumping. If this pumping is replacing 
atoms at the same rate as they are being removed then it may be possible for the rate of evolution of the spatial mode 
function of the condensate atoms to be much slower than the rate of coupling. In this case the above approximation 
will become more accurate. 

In this paper we proceed by first assuming that the mode function and the oscillation frequency are constant and 
then determine the memory time. This memory time is then compared to the time scale of the dynamics of the 
trapped mode of interest. If the memory time is much shorter than the time scale of the mode dynamics the above 
assumptions are valid. If this is not the case then we have a fully non-Markovian decay with a time-varying system 
frequency and coupling constant. Obviously it is of interest to determine the regions of validity of the two cases. 



C. Green's Functions 



To determine the validity of the Markov approximation we need to calculate fm{T~) in the presence of the potential 
V{j (x) and for an interaction region k(x) . A convenient way to do this is via the Green's function introduced above. 
The single particle Green's function for the untrapped particles, Eq. (|l3|), can be written in terms of path integrals 
(see for example, pif ) as, 

G(x, t; x', f ) = £* dx(r) exp j X (18) 

where t = t — t' and where S is the action of the particle given by, 5[x(r)] = f., drL(x, and L = ^m(^) 2 — 
V£ s (x.) is the La grangian for the untrapped particles. We are interested in the case where the Lagrangian is the 
sum of the Lagrangians in each dimension. In this case the Green's function factorizes into three one-dimensional 
Greens functions. In general, the path integral is difficult to calculate. However, one can make a semi-classical 
approximation to the Green's function. It turns out that this approximation is exact for potentials up to quadratic 
order in the coordinates. In fact, the semi-classical approximation is justified by approximating a particular potential 
by a quadratic potential f|24j|. In this work we will only deal with quadratic potentials. 



In the following sections (Sec. Ill and IV) we will consider interaction regions that are independently Gaussian 



shaped in all three dimensions, so that k(x) = k(x)k(j/)k(z), where each k(J), for j — {x, y, z}, takes the form 

«(i) = -^-rexpJ-i^l, (19) 

where Oj is the width of the Gaussian in the jth dimension. Note that we are ignoring any multiplicative constant of 
k(x) as we are interested in the relative fall off of the memory function. The assumption of a Gaussian interaction 
region allows us to calculate memory functions exactly in many situations. We do not expect the exact shape of the 
interaction region to affect the order-of-magnitude estimates for the memory time that we make in this paper. 



D. Properties of the output-coupled atoms 

In the previous sections we have made some general considerations about the Markov approximation for the damping 
of the trapped mode. There is, however, another aspect to the problem; that of determining the properties of the 
output-coupled atoms. For example, let us assume that to a good approximation the system exhibits Markovian 
damping and the atoms leave the interaction region with a reasonably well-defined momentum. Due to the dispersive 
nature of the vacuum for atoms the properties of the untrapped field will have a non-trivial dependence on position. 
Very close to the trap, atoms will not have traveled very far and dispersive effects may not be large. However, if 
the atoms experience a lot of dispersion then the properties of the field will correspond to properties of the trapped 
mode averaged over some time. Let us assume that there is a measurement device localized about a position xo that 
is making destructive measurements (destructive in the sense that the detector scatters atoms into free modes far 
from those of interest; examples of such detectors are a hot wire or ionization by a laser) on any atoms that interact 
with it. In this case there will be an uncertainty in the time of emission of an atom that is detected at Xo- The 
description of a continuous measurement process becomes much more complicated in this regime p5| ]. We call this 
the non-Markovian regime for the measurements. 

To describe such a situation (see Appendix A) we can define a response function of the system to a particular 
measurement device, (in analogy with the memory function) , as 
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h x (t-t')= / dxdx / x(x-xo)K(x',i / ) G ( x ^; x ': i ')> 



(20) 



where %(x — x ) describes the spatial extent of the detector. Ideally, the spatial extent of the measuring device will 
be smaller than that of the interaction region, otherwise much of the uncertainty will be introduced by the detector 
itself. h x (t —t') is the probability amplitude for a particle that is emitted in the interaction region at time t' to be 
detected at time t by the detector. 

We can define a memory time for the detection as the time interval between the earliest and the latest time that 
a particle could have been emitted. The memory time for measurements of the output can be analyzed in the same 
way as for the damping. If the memory time corresponding to this response function is much shorter than the time 
scale of the system dynamics then we can make the Markov approximation for the measurements. If this holds then 
a detection time can be considered to correspond exactly to an emission time and all the moments of the measured 
field are proportional to those of the trapped mode. This will often be a stronger condition than that for Markovian 
damping. 



III. RADIO-FREQUENCY OUTPUT COUPLER 



In the radio-frequency output coupler a radio wave of frequency uj r f induces transitions between trapped and 
untrapped (or anti-trapped) magnetic sub-levels of the atoms. The strength of the coupling is given by the Rabi 
frequency, f2 = #£tBohr | -B | / \/2h, written here in terms of the magnetic field B and the Lande g-factor. The waist of 
the r.f. wave is assumed to be much broader than the spatial mode function of the trapped mode and so from Eq. (|^) 
the interaction region becomes re(x, t) = u a (x)e wt , where, if a is the condensate mode, u a {~x) is determined by the 
solution to the time-independent GPE. However, we assume here that it is valid for our purposes to approximate this 
mode function by a Gaussian. The energy difference between the untrapped level and the center of the trap is given 
by Vq — Vt(0). The untrapped atoms are free to propagate away and are in general subject to accelerating potentials, 
V^ ff (x). The general situation is depicted in Fig.|l|. 

In the following sections we will determine memory times for a number of relevant potentials for the untrapped 
atoms. 



A. Free space 



In order to emphasize the effect of the external potentials we first consider the case when the atoms are coupled 



into free space. In the x dimension the free space, (V^ ff (x) = 0), Green's function is 



G{x,x';r) = _l=exp(*^-^ \ , (21) 



where A = h/2m is a measure of the rate of spreading of the wave packet. Integrating this over the Gaussian 
integration region and multiplying the three integrals for each dimension together yields the memory function 

f m (r)=Y[A 3 (r)e-^, (22) 

3 

where j = {x, y, z}, ojq = /.t + Vq — cj r f and, 

A,(r) = g ■ (23) 

V°i + iXt 

A Gaussian wave packet will keep its Gaussian shape (in real space) but will increase in width over time due to the 
fact that it contains a range of velocity components. The overlap of the wave-packet with itself as a function of time 
is given by the memory function. 

The radio frequency field couples an atom from the trapped mode to modes of the untrapped field with frequencies 
around u>q = Vq + /i — w r f so that v = ^j2?iujQ/m is the magnitude of the mean velocity of the output coupled atoms. 
When the radio frequency field is on resonance with the trapped mode, (luq — 0), particles are coupled out with a zero 
mean velocity. The particles can only leave the interaction region by quantum mechanical spreading of their wave 
packets. The memory function will then decay as t~ 3 / 2 for long times, r ^> cr|/A. If the output coupled particle has 
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an initial mean velocity, wo ^ 0, the memory function will decay at a faster rate, as the atom will leave the interaction 
region more quickly. The case where the mean kinetic energy of the atom is much higher than the coupling rate, 
u>o ^> 7, is very similar to the optical case and a memory time can be defined for which 7 _1 ^ T m 3> luq 1 - 

Let us investigate the memory time more quantitatively. In order to compare memory functions which decay in 
very different ways we consider the ratio of the magnitude of the two integrals in Eq. ( M ) , 



R = 



t-T m 

dsf m (t - s) 



t 

dsf m (t - s) 



(24) 



The quantity R is a measure of the inaccuracy of the approximation that a particle has left the interaction region for 
a particular choice of T m . Often we are interested in the inverse, i.e., the memory time given a certain lower bound on 
the accuracy of the approximation. Doing the integrals in Eq. (^) we can find the ratio, R, in terms of the memory 
time (see the Appendix B for details of the calculations). First consider the case when the system is on resonance, 
luq = 0. In the symmetric interaction region case, a = <jj, the ratio reduces to 

R= K + (AT m )2]i/4- ( 25 > 

Inverting this equation we find the memory time in terms of the ratio, T m > a 2 / ' \R 2 , where we have assumed R <C 1. 
In the asymmetric case, the memory time is the same as the symmetric case where a is the broadest width of the 
interaction region. This memory time depends on the square of the width of the Gaussian and so is a sensitive 
function of the size of the interaction region. In experiments performed to date the size of a condensate in the 
broadest dimension has been a > 10/im. For 87 Rb [m ~ 10~ 25 kg), if we assume a ratio of R = 1CP 2 and a size of 
a ~ 10/im, this already gives a very long memory time of the order of 10 3 s. 
When the atoms are coupled out with an initial velocity, ujq ^ 0, we get 

^ > -^r, (26) 

\Rd 

where d is the dimension of the untrapped field, such that, in the case where the interaction region is cigar shaped, 
(e.g., a — a x — a y and a -C cr z ), d = 2, and when it is pancake shaped, (e.g. a x = a y , and a x ^> a z = a), d = 1. 
The initial velocity produces an oscillation of the memory function which when averaged over many oscillations leads 
to a reduction in the memory time compared to the on resonance case. For a cigar shaped region with a ~ 10/im 
and R — 10 -2 this gives a much shortened (but still relatively long) memory time of T m ~ 10s for Rb atoms. These 
calculations show that atoms of low velocity leaving the interaction region by the spreading of their wave packets 
linger in the region of interaction for times much longer than Is. 

For very weak coupling, the oscillation frequency, u>o, due to the initial kinetic energy of the atoms may be much 
greater than the coupling rate 7. The memory function averaged over times much longer than the time scale defined 
by the damping 1/lu c <C 1/7 is given by 

sinh([cr 2 + iAr] !f) 



/m(r) « VL „ - J x > . (27) 

<j z + i\t 

An integral over this function converges, as it acts like a sinc(w c r) function for large r 3> a 1 j\ and we can define 
a memory time by T m ~ 1/lo c . Note that we do not analyze the form of the decay in this case as it is simply due 
to our choice of a sharp cutoff to restrict the frequencies (see Appendix B). In this regime, there is essentially no 
difference between the Markov approximation for the damping in the optical and the atomic case. This is due to the 
fact that we have implicitly assumed that the coupling constant is approximately constant across the frequencies of 
interest and that we are on a linear part of the dispersion curve or, equivalently, that the atoms have a limited range 
of velocities about a fast mean velocity and so have fast propagation times across the region of interaction. 



B. Gravity 

In most situations atoms coupled out of a trap will be subject to gravitational forces. It is therefore of interest to 
consider the effects of gravity on the length of the memory time. 
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The Green's function in the z-dimension for a gravitational potential (z) = —gmz has the form 



ni ' \ I 1 \ \ z - z 'f . gm{z + z') mg 2 T 3 \ 

The first term in the exponential being the usual dispersion term and the second and third terms can be recognized 
as phase shifts due to the potential and kinetic energies, respectively. Integrating over a Gaussian shaped interaction 
region and assuming free space Green's functions for the other two dimensions the memory function becomes 

Ur) = IlA^exp {- p?f) 2 - i^r* - i. r\ , (29) 

3 

where Tlluq = Tifi + Vq — hui r f is the initial energy of the output coupled particles (measured from V^(z = 0) =0). 
Under the influence of gravity an initial Gaussian wave packet preserves its Gaussian shape (in real space) but the 
peak of the Gaussian propagates at a velocity v = gt after a time t in the — z direction. The Gaussian decay of 
the memory function (given by the first term in the exponential) is due to the gravitational potential accelerating 
particles out of the interaction region. If we assume that this Gaussian decay is the dominant process for short times 
the memory function can be written as 

/ m (r) ocexpj-^j, (30) 
where a T = \/2h/mga z . The ratio, R, defined by Eq. ([hi]), becomes 

(3D 

V7rT m 

where we have used the asymptotic behavior of erfc(a;) [p6[ . Inverting this, and assuming T m 3> er r , we find T m > 
(7 T -\/21n(l/i?) = 27i-y/ln(l / R) / mga z . In this case, T m is inversely proportional to the size of the interaction region 
in the z-direction. Surprisingly, the time for a particle to leave the interaction region actually gets shorter as the 
interaction region gets larger. For an interaction region of size a z ~ 10/iin this Gaussian envelope gives a memory 
time of T m ~ 10 _5 s for Rb and T m ~ 10~ 4 s for 23 Na and will becom e shorter for a larger interaction region. The 
memory time is largely independent of R for R -C 1, as, in this case, \/\n(l/R) ~ 1. 

Let us now consider the second term in the exponential in Eq. (|29|). This term describes the property that after 
a certain time the particle is traveling at a velocity high enough that it causes a rapid oscillation of the memory 
function. This oscillation will average to zero over a time scale much longer than the oscillation frequency. It makes 
physical sense to estimate a memory time as the time it takes a particle, accelerating from rest under gravity, to 
reach a velocity that is high enough that the particle can no longer be coupled back into the trap, t = v/g, where 



v ^S> y/2hr)jra. From Eq. (29) the memory function will be oscillating much faster than 7 if T m is such that 



2^»7- (32) 



Rearranging this, we find T m 3> 2/gy/jh/m which supports our initial estimate quite well. Note that this memory 
time is independent of the size of the interaction region. For Rb and assuming the oscillation is 10 times greater than 
7, we can estimate the memory function in terms of 7 as T m ~ 10 -5 ^/t~s 3 / 2 . 

The particle will leave the interaction region sooner than it can gain the required kinetic energy if a T is less than 
the time it takes to make one oscillation a T < (iSn/mg 2 ) 1 ^ 3 . We can rewrite this condition as a z > {A&ir^gi:)^ 1 / 3 . 
For Na, a z would need to be greater than 6/im and for Rb, a z needs to be greater than 0.1/im for the oscillation 
to become important. These values correspond to very small condensates and therefore we expect the memory time 
to be determined by the time it takes the particles to leave the interaction region in most situations. A plot of two 
possible situations is given in Fig.||. 

C. Collisional repulsion and the anti-trapped case 

If particles are being coupled out of the condensate mode into an untrapped state they will see a repulsive potential 
due to the condensate atoms left in the trap (in the repulsive interaction case) and the shape of this potential will be 
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proportional to the density distribution of the condensate. In the Thomas-Fermi approximation the condensate density 
takes the shape of the trap potential f27| , UN\u a (x-)\ 2 = max [// + Vq — I^r(x), 0]. In the present experimental situa- 
tions the traps have been harmonic in all three dimensions. The repulsive potential for the untrapped atoms will then 
be a three-dimensional inverted harmonic potential, given by V^ ff (x) = eUN\u a (x.)\ 2 — max[e/i — \vn ^2 e(wj) 2 j 2 , 0], 
where e is the ratio of the scattering lengths between a trapped-untrapped atomic collision and a trapped-trapped 
collision. In this paper we assume that in the region of interaction we can approximate the effect of the cutoff inverted 
harmonic potential by a inverted harmonic potential that is not cutoff and is therefore quadratic everywhere. In 
general, the Green's function depends on the potential everywhere, but we are only interested in the time until a 
particle is repelled from the interaction region and as long as the particle energies are not too close to the cutoff we 
can neglect the global effects due to the shape of the potential outside the interaction region. 

In addition, if we are considering particles that are output coupled into spin states that are repelled by the trap this 
is also an inverted harmonic potential. In this case the repulsive potential will be V^ ff (x) = eUN\u a (x)\ 2 — Vr(x), 
where we have also included the collisional repulsion. In this case, e is the ratio between a trapped-antitrapped atomic 
collision and a trapped-trapped collision. 

The Green's function for an inverted harmonic potential is easily determined from the Green's function for an 
harmonic potential p4| with the substitution ui — > iui. Along a single axis of the inverted harmonic potential we have 



G(x, x'\ t) — 



mui x 



2nih sinh w t t 



exp ■ 



imuo x 



\(x 2 + x' 2 ) cosher — 2xx'] — i-^-r \ , 
2hsmh uj x t n \ 



(33) 



where Vq is the potential at the center of the inverted harmonic. This is the Green's function of the untrapped atoms 
if we can assume that the effect of the inverted harmonic potential is much greater than that of gravity. In light of 
the above discussion it is rather inconsistent now to assume a Gaussian interaction region. However, we do not expect 
the exact shape of the interaction region to effect our results dramatically. This is borne out by results which we will 
present below. 

The corresponding memory function in the Gaussian case becomes 



/m(r)=n 



coshc 



A 



A 



sinh Ljj t 



-1/2 



(34) 



where luq = fi + Vo/h 



uj t { — Vq/Tl. For times t ^S> this becomes 

1 -1/2 



/ m (r) c J; 



1+i 



A 



A 



cxp j- yT - icj r| 



(35) 



This exponential decay describes particles being repelled out of the interaction region by the inverted harmonic 
potential. Assuming this exponential decay is the dominant process for short times we can estimate the usual ratio of 
integrals by R ~ exp(— 3wT m /2), where to = (lo x +u y + u> z )/3 is the mean of the inverted harmonic trap frequencies. 
Inverting this we get T m > (2/3w) ln(l/i?), which is independent of the size of the interaction region. For very 
asymmetric traps the trap frequency with the largest magnitude will define the memory time as particles will be 
repelled fastest in this dimension. 

In the above calculations we assumed a Gaussian interaction region. However, these results hold for general 
interaction regions. This can be seen by taking the limit r ^> 1/oJj (as suggested by Eq. (|3^)) of the Green's function 
itself, 



G(x,x';t > l/ujj) 



moj x 
2irih 



exp 



im,LO 



27) 



-(x 2 +x' 2 ) 




(36) 



where we have assumed 1 3> e _Wxr . This decays as an exponential with increasing time, independent of the shape of 
the interaction region. 

In the untrapped case where the potential is determined by collisional repulsion ojj = \/tuJ and Vq — efi. Assuming 
that e ~ 1 we can estimate the memory time (as defined by the exponential decay) as simply one over the mean of the 
trap frequencies. Typical trap frequencies range from 2it x 10 — 2n x 400 Hz. From the experimental parameters of || 
we can determine the value uj = 2ir x 219Hz for the average trap frequency which gives a memory time of T m ~ 1ms. 

In the anti-trapped case we have uij — (\/e+ an d Vq = efj, — Vq. The effective inverted harmonic trap potential 
is estimated to be approximately twice the strength of the trap potential as it is the sum of the repulsive potential and 
the collisional repulsion. So the memory time calculated for the untrapped case will be halved due to the repulsive 
trap potential. 
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D. Summary of the radio- frequency output coupler 

A summary of the calculated memory times for the radio-frequency output coupler is given in Table ffl, where R 3> 1 
and d is the dimension of the space of the untrapped atoms, e.g., d = 1 corresponds to a cigar shaped interaction 
region. 

If particles are coupled out with a large kinetic energy, loq 3> 7, then a memory time can be defined by 7 -1 3> 
T m 3> U!q , in the same way as for an optical system. In the optical case the field inside a cavity is coupled directly 
to the field outside and particles tunneling out of the cavity have nearly the same energy as they had inside the 
cavity. However, in the atomic case the coupling between the trapped field and the untrapped field is mediated by 
the radio-frequency field. Particles that make the transition between a trapped state and an untrapped state have 
an initial energy equal to the detuning between the frequency of the r.f. field and the energy difference between the 
two states, see Fig.|l|. This means that in the case where ojo is large the radio-frequency is far from resonance with 
the condensate mode. If a populated non-condensate mode is close to resonance the particles in the output will be 
thermal particles. In order to maximize the coupling to the condensate, the r.f. field will need to be on resonance 
with the condensate mode and consequently we expect particles to come out with a small mean energy, luo w 0. The 
slow moving particles are then subject to the accelerating potentials of gravity and collisional repulsion. 

From our calculations, gravity seems to be the force that determines the memory time. It yields a memory time of 
T m ~ 10 -2 ms for a Rb condensate with a width in the z-direction of the order of 10/im and decreasing with increasing 
width. Sodium, being lighter, has a longer memory time. Whereas collisional repulsion gives a memory time of the 
order of the inverse trap frequencies T m ~ 10 _1 ms. This depends on the scattering length of a trapped-untrapped 
collision and is independent of the size of the region of interaction. 



IV. RAMAN OUTPUT COUPLER 

We can model the case of a Raman output-coupler by replacing the coupling constant in the direction of the kick, 
k(x), by 

k(x) -» K(x)e- lAx , (37) 

where HA is given by the difference between the momentum of the two photons involved in the Raman transition, 
A = fci - k 2 . 

Assuming the Raman kick is the dominant process, the Green's function is the free space Green's function. The 
memory function, assuming a Gaussian shaped interaction region, is then 



/ m (r)=nA,(r)ex P {-l^^-^} 



(38) 



where in this case loq = ji — Aw -\-Vo/Ti, and we have defined Aoj = lo\ — L02 as the difference between the frequencies 
of the two lasers. 

In this case it is possible for focused laser beams to be of the order of the size of the condensate. The width of 
these beams would then define the size of the interaction region. However, focusing the beams to a small region will 
tend to lead to diffraction. The momentum of the kick will not be well defined in this case. A comparison between 
the situation when the momentum kick is well defined and when it is not is shown in Fig.||. If the momentum kick is 
not very large compared to the characteristic inverse length of the interaction region, A ~ l/cr x , then for long times 
the behavior is similar to the radio frequency case (i.e., repulsive potentials will dominate). However, where the kick 
is much greater than the inverse size of the region A 1 j a x we can expand the term in the exponential in Eq. (J3S~ 
and write the memory function as 



/ m (r) » n A,(0) exp {- ^ - z(AA^ + w )r| , (39) 
j T 

where <j T — er x /\/2AA. This will dominate for time scales r <C o^/A. In the case when the kick is well defined 
the memory time has a very simple classical interpretation in terms of the time a particle takes, t = d/v, to cross 
a distance d = 2a going at velocity v = HA/m. The first term in Eq. ([59]) is a Gaussian envelope of width a T and 
represents the overlap between the wave packet of a particle propagating away and the interaction region. As in 
the gravitational case this Gaussian envelope determines a memory time in terms of the usual ratio of integrals as 
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T m > a x y/ln(l/R) / AA . In the best-case scenario, if the two laser beams involved in the Raman transition are 
counter-propagating, then at optical frequencies the kick given to the atom will be of order A ~ 10 7 m _1 . This yields 
a memory time of T m ~ 10~ 3 s for interaction regions of a x = 10/im for Rb atoms. 

The second term in Eq. ( |39|) is recognized as a kinetic energy term for a particle propagating with momentum 
h(A 2 /\ + ujq). If the kick is large the particle will be given a large kinetic energy and a memory time can be 
determined by a high oscillation frequency. If A 2 A + ujo ^> 7, then we can determine a memory time as T m ~ 1/cj C) 
where A 2 A + u>q ^> u> c 3> 7. Unlike the radio frequency output coupler the Raman output coupler (in the case where 
the kick is well defined) is very similar to the laser in that condensate atoms will leave the trap with a finite, and 
possibly large, kinetic energy in a well-defined direction. If the kick is large enough the effects of gravity may be 
negligible for short distances. 

DISCUSSION 

We have investigated the regimes of validity of the Markov approximation for atoms that are being coupled out of 
an atomic trap by an output coupler by determining the memory times (or correlation times) of the output coupled 
atoms. A memory time can be defined by the time taken for a particle to leave the region of phase space of the 
untrapped field where it can be coupled back into the trap. After this time, if a particle has not been coupled back 
it can be considered to have irretrievably left the trap. 

Atoms coupled out by an r.f. field that leave the interaction region by the free space spreading of their atomic wave 
packet alone can have long memory times, (T m 3> Is), which are dependent on the square of the size of the broadest 
dimension of the interaction region. These memory times are reduced if atoms are coupled out with a large mean 
kinetic energy. However, in this case the coupling will be far from resonance with the condensate mode. For the r.f. 
output coupler, gravity (which will nearly always be present) dominates in many cases, and yields a memory time 
that ranges from 10~ 2 — 10 ms for typical experimental parameters. Collisional repulsion leads to a memory time 
of the order of the inverse mean trap frequency, Q, and depending on the ratio of the scattering lengths between a 
trapped-untrapped and a trapped-trapped collision, e: T m ~ which is slightly longer than the memory time 

for gravity for typical experimental parameters. 

Raman output coupling has a memory time which depends on the inverse strength of the momentum kick, HA, 
given by the light to the atoms and on the size of the condensate in the direction of the kick, a. The kick must 
be well-defined, A ^> 1/cr, to produce a reduction of the memory time. If the light beams are focused too tightly 
(A ~ l/c) diffraction effects will dominate and the Raman output coupler, like the radio-frequency output coupler, 
must rely on external potentials such as gravity to determine a memory time. In the best case scenario, two counter 
propagating laser beams will produce a kick of A ~ 10 7 m -1 giving a memory time of T m ~ 1ms. The Raman output 
coupler has the nice property that in the regime where a short memory time is produced the atoms are given a 
well-defined momentum kick, producing a beam of atoms [Q . 

These memory times must be short compared with the time-scales of the system evolution, and the coupling 
rate 7 -1 , in order to make the Markov approximation. The most important upper time limit on 7 -1 is the correlation 
time of the condensate, r c . 7 -1 must be short compared to r c otherwise the output-coupled atoms will not be 
correlated with each other. Current experimental estimates for the correlation time give a range of from 0.1 — Is, i.e., 
long compared to the calculated memory time corresponding to gravity. Therefore, the Markov approximation will 
be valid for a range of coupling rates where r c 3> 7 _1 ^> T m . 

When the Markov approximation is valid a Markov master equation may be used to solve for the evolution of the 
trap modes undergoing damping. When the Markov approximation cannot be made new methods must be employed 
f3|, P5|]2gf| , and novel behavior will be observed . 
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APPENDIX A 



A localized system (in this case the trapped atoms) interacting with a bath (the modes of the untrapped field) can 
be thought of in terms of inputs and outputs to the system [ p9| . The initial bath propagates towards the system; 
interacts with the system; propagates away again; and is eventually measured. To formalize this idea, input and 
output fields are defined by 

iMx, t) = Ul(t, io)V-c/(x, t )Uu(t, t ), (40) 

where to < t and is usually taken to be in the distant past, and 

</w(x,<) = Uu{h,t)ipu(x,h)Ul(ti,t), (41) 

where t\ > t and is usually taken to be in the distant future, and where the field operators are in the Heisenberg 
picture. Ujj is the evolution operator defined in terms of Hu alone. The relation between the input and output fields 
and the trap mode is given by 

4ut(x,t) =^ in (x,t)+V7 / ds[^ in (t,x.),i\s)}a(s). (42) 

Jt 

The output has a contribution from the input fi eld a nd the trapped mode at earlier times. 

We assume that our measuring device (of Sec. II D) measures normally ordered moments [p0| of the quantity 

*out(*) = ydxx(x-x )Vw(x,i), (43) 

where x( x — xo) describes the spatial extent of the detector. We are assuming that the actual measurements take 
place over a very short time. The response function for the system given a particular measurement device is defined 
as 

h x (t-t') = / dx X (x - x )[4.(i, x), £ f (*')]■ (44) 



h x (t — t') is the probability amplitude for a particle that is emitted in the interaction region at time t' to be detected 
at time t by the detector. This becomes more obvious if we write it in terms of the Green's function for the untrapped 
field Eq. @. 

We can define a memory time, T m , as the time interval between the earliest time a detected particle could have 
been emitted and the latest time the particle could have been emitted. T m exists if we can write 

t pt 
dsh(t — s)a(s) ~ / dsh x (t — s)a(s), (45) 



t-Tm 



where we have ignored any constant time delay that the response introduces as it will not effect steady-state results. 
If this memory time is much shorter than the time scales of the system dynamics we can write 



# out (t) ~ % a {t) + VW). (46) 

where 7 is defined in a similar way to 7' for the damping and ^i n (t) is the contribution of the input field to the 
output. If this holds a detection time corresponds exactly to an emission time and all the moments of the measured 
field are proportional to those of the trapped mode. 



APPENDIX B 



The ratio of the magnitude of the integrals over the memory function from t — T m to —00 and from t to —00 in the 
free space case is 



R = 



t-T m 



dsJ|Aj(t - s) 



/ ds JjAj(i - s)e" 
J — 00 



(47) 



■iu)Q (t—s) 
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Doing the integrals in this equation allow us to determine the ratio in terms of the memory function. 
In the symmetric interaction region case, where a = cr 3 ■ , we can do the integrals in Eq. (|47| ) and we find 



R 



ere 



Vo 3 



i\T„ 



A 



"erfc 



LU 



A 



A erfc 



(48) 



where erfc(z) = 1 — erf (2) is the complementary error function. In the case where ujq = this ratio reduces to Eq. (|2E 
In the asymmetric case with ujq = 0, although the intermediate behavior is governed by the narrowest dimensions, the 
memory time is determined by the broadest dimension. This is due to the fact that in strictly one and two dimensions 
we cannot define a memory time as in these cases the integrals in Eq. ( pi] ) diverge. This is a fundamental property 
that arises from the 1/y/r dependence of the one-dimensional Green's function. In conclusion, the third dimension is 
required for the system to be dissipative in the luq — case. 

On the other hand, the case ojq ^ is important in the asymmetric case as it is well known that a multiplicative 
oscillating factor can make an otherwise divergent integral convergent. In the case when w / we can use the 
asymptotic expansion for the error function, [Eq], y/nze z erfcz ~ 1 + l/2z 2 as z — > 00 for |arg(z)| < 37r/4, to show 
that for long times T m > l/w , R oc 1/[ct 4 + (AT m ) 2 ] 3 / 4 . For u > A/ct 2 , Eq. (El simplifies to 



R 



]CT 4 + (AT m ) 2 ] 3/4 ' 

In the case where the interaction region is cigar shaped, e.g., a = a x = a y and a <C ct 2 , the ratio becomes 



(49) 



R = 



{c (4+zT m )} 



(50) 



where Ex (z) is the 1st order exponential integral 



Ei{^°-} 

I . Note that a one-dimensional interaction region corresponds to 



a two dimensional bath and vice versa. For long times we can write R oc a 2 / \J cr A + (AT m ) 2 , where we have used the 
asymptotic expansion Ei(z) » e~ z /z as z — > 00 for |argz| < 3ir/2. The ratio simplifies in the ujq ^> A/ct 2 limit to, 



R 



+ (AT m ) 2 ]V2' 



(51) 



for all times, T m . The memory function is then given by T m > a 2 /RX. 

If the interaction region is pancake shaped, e.g., a = a x , a y 3> ct, and ct 2 3> ct we get 



R = 



erfc < 


V"o(T+ ir -)} 


erfc < 




\ 



(52) 



When ujq A/ct 2 this becomes 



R 



K + (AT m ) 2 ]V4- 



(53) 



Yielding a memory time with the same dependency as the low frequency symmetric case. Putting this all together 
yields Eq.©. 

If the atom is coupled out with a large kinetic energy it is useful to consider the memory-function written in the 
form of an integral over frequency, 

/•oc 

/m(t-0=/ duD^K^e-^-W-^, (54) 
Jo 

where to — u>k, D(ui) is the density of states and k(oj) is found by transforming the effective coupling constant to 
k-space (where k is the label of the modes of the untrapped field) and then using the dispersion relation to make a 
change of variables to frequency space. 
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In dealing with averages over oscillatory functions it is necessary to specify a time scale over which averages are to 
be taken. In frequency space this corresponds to considering only a range of frequencies in the integral in Eq. (^4|). A 
physical frequency-dependent coupling, k(lo), will naturally limit the range of frequencies and we can consider memory 
times defined by Eq. (jlj) . 

The time scale that we are ultimately interested in though is 7 -1 . We should therefore average over time scales 
short compared to 7 -1 . In frequency space, we can introduce a simple cutoff frequency uj c 3> 7 so that the function 
is averaged over times of order w" 1 , 



f m (t - t') = e-^ 1 / " duD(u + u )\K(u + w )\ 2 e- iu <- t - t '\ 



(55) 



where too = /u + v > lu c . Obviously, loq must be much greater than this for the function to be attenuated with 
increasing time. We are essentially band filtering the memory function so we refer to this as the filtered memory 
function. Physically, we are neglecting particles with energy greater than h(u>o +to c ) and less than h(coo ~oj c ) because 
they are a long way from resonance. 

We can define a memory time in terms of the filtered memory function f m {t — t') by 



dsf m (t 



t-T„ 



> 



t-Tm 



dsf m (t 



(56) 



In many cases T m ~ l/u> c - Note that a memory time defined in this way relies on the fact that we have assumed 
^0 > 7- 

The memory function averaged over times l/w c where lu 3> lo c 3> 7 is given by Eq. (|55|). In the symmetric case, 
D(uj) oc ^/To and for a Gaussian shaped interaction region, |k(w)| 2 = exp(— ^-w). The filtered memory function 
simplifies for long times to 



/ m (r) oc sinh (Jcr 2 + iXr] 



1 



a 2 + iXr 2{a 2 + i\r) 2 



(57) 



where we have again used the asymptotic expansion for the error function. The second term will tend more rapidly to 
zero and so the first term will define the memory time. An integral over this term converges as it acts like a sinc(w c r) 
function for large r 3> a 2 /A and we can define a memory time by T m ~ l/w c . Note that we do not analyze the form 
of the decay in this case as it is simply due to our choice of a sharp cutoff to restrict the frequencies. 

If the interaction region is effectively two dimensional then the density of states becomes D(lj) oc 1/s/uj with |k(o;)| 2 
the same as above. The filtered memory function for long times also simplifies to 



/m(r) oc 



sinh ( [a 2 + iXr] 
a 2 + iXr 



(58) 



where we have again used the asymptotic properties of the error function. In a cigar shaped interaction region the 
density of states is D(u>) oc 1 and the filtered memory function has the same long term behavior as the two and three 
dimensional cases. 
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x (ms) 

FIG. 2. Plot of the real part of the memory function for a particle in a gravitational potential. The memory function is 
plotted for Sodium atoms with a symmetric interaction region a x — a y = a z and the frequency loq = 2tv x 100Hz. The solid 
line corresponds to a z = lO^im where the time for a particle to leave the spatial interaction region determines the memory 
time. The dashed line depicts the case of a smaller interaction region a z — 1/im where it is possible that the oscillation due to 
the accelerating particles velocity determines the memory time. 
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TABLE I. Summary of memory times for the radio-frequency output coupler. 
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FIG. 3. This figure is a comparison between the case when a particle is given a well defined momentum kick and when 
it is not for the Raman output coupler. We have plotted the real part of the memory function as a function of time for the 
parameters A = 10 6 m -1 , a y = a z — 10/xm and loo = 2ty x 100Hz. The solid line is the case a x = 10/im so that A ^> l/a x . The 
dotted line is the case a x = 1/J.m, so A ~ l/a x . In this last case the memory function does not exhibit the Gaussian decay as 
the momentum kick is not well defined. 
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